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ABSTRACT 

Let  the  points 

(1)  (i  *  1, . . . ,kj  k  >  2),  a  <  x1  <  x2  <  •••  <  xk  <  b  , 

I  *  [a,b] ,  (-•  <  a  <  b  <  •) 

be  prescribed.  Furthermore,  let  m  and  n  be  integers  such  that 

1  <n<k<m, 

a  m 

and  define  the  polynomial  class 

II  =  (p<x);  P(x)  e  »  ,  P(xH)  *  y,  ,  (i  =  1,...,k)}  . 
w  nil 

Within  II  we  determine  P_(x)  as  the  solution  of  the  extremum  problem 

m  in  *■ 

(2)  /  (P(n,(x))2dx  -  minimum  for  P(x)  e  Mp,  . 

I 

Finally,  let  S(x)  «  S2m_^(x)  be  the  natural  spline  interpolant  of  degree 
2n  -  1  of  the  k  points  ( 1 ) .  Our  main  result  is 

Theorem- 1.  1.  There  is  a  unique  polynomial  Pffl( x)  which  is  the 

solution  of  the  minimum  problem  (2). 

2.  We  have 

lim  Pm(x)  =  S(x)  uniformly  in  x  e  I  . 
m+0° 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  finite  interval  [a,bl  we  have  prescribed  abscissae 
<  *2  <  •••  <  and  corresponding  ordinates  y^,y2»  Let 

Six)  =  S2k_1(x)  be  the  natural  spline  of  degree  2n  -  1  that  interpolates 
chose  k  points.  This  requires  that  1  <  n  <  k.  Furthermore,  let  m  be  an 
integer  such  that  m  >  k,  and  let  Pm(x)  be  the  polynomial  of  degree  at 
most  m  that  interpolates  the  k  points,  and  such  as  to  minimize  the 
integral 


/  (P<n)(x))2dx  , 
a 

within  the  entire  class  of  polynomials  of  degree  m 
points.  It  is  shown  that  as  m  ♦  the  polynomial 
spline  S  ( x)  . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


INTERPOLATING  SPLINES  AS  LIMITS  OP  POLYNOMIALS 
I.  J.  Schoenberg 

For  Alexander  Ostrovski  on  his  90th  birthday  on  September  25,  1983, 
from  one  of  his  grateful  students. 

1 .  Introduction.  Let  the  points 

(1.1)  (xi,yi),  (i  -  1,...,kj  k  >  2),  a  <  x1  <  Xj  <  •••  <  x^  <  b  , 

I  -  [a,b] ,  (-•  <  a  <  b  <  «)  , 

be  prescribed.  The  basic  interpolant  is  the  Lagrange  interpolating 
polynomial.  If  additional  consecutive  derivatives  at  the  points  (1.1)  are 
available,  we  can  construct  the  Hermite  interpolation  polynomial.  In  the 
absence  of  such  additional  data,  we  propose  here  the  following  construction: 
Let  m  and  n  be  integers  such  that 

(1.2)  1  <  n  <  k  <  m  , 
and  let  us  consider  the  polynomial  class 

(1.3)  II  -  {P(x)f  P(x)  e  *  ,  P(x  )  -  y .  ,  (i  -  1 , . . .  ,k)  }  . 

m  rail 

Within  this  class  we  determine  the  polynomial  Pm(x)  which  is  "most  nearly  a 
polynomial  of  degree  at  most  n  -  1  in  the  interval  I."  We  interpret  this 
requirement  to  mean  the  Pm( x)  is  the  solution  of  the  extremum  problem 

(1.4)  /  (P*n*(x))2dx  -  minimum  for  P(x)  e  II  • 

I  m 

Equivalently:  Writing 

(1.5)  -  inf  /  (P(n,(x))2dx  , 

'  pen  i 
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the  polynomial  Pm(x)  ia  uniquely  defined  by 


(1.6) 


/  (Pin)(x))2dx  -  Pn(x)  e  1^ 


a 


Of  course,  the  existence  and  uniqueness  of  Pm(x)  is  yet  to  be  established. 

Our  main  subject  is  the  behavior  of  ?m(x)  as  m  +  ».  The  statement  of 
our  result  (Theorem  1  below)  requires  some  known  properties  of  natural  spline 
interpolation.  We  describe  its  definition  and  the  three  properties  that  we 
need. 

I.  Let  n  be  an  integer  such  that 

(1.7)  1  <  n  <  k  , 

and  let  S(x)  «  s2n-1*x*  denote  a  function  satisfying  the  following  four 
conditions  t 

1»  s(x)  e  c2n"2(»)  , 

2*  S( x)  e  ir  in  each  interval  (Xj,x^+1),  (j  -  1,...,k  -  1), 

3#  S(x)  e  it  «  in  (-*,x_),  and  S(x)  e  it  in  (*,+»)  , 
n- 1  —  1  '  n- 1  —  x 

4*  S(xi)  -  yif  (i  »  1,...,k)  . 

Then  S(x)  ia  uniquely  defined  by  the  conditions  1*  to  4*. 

The  function  S(x)  ia  called  the  natural  Bpllne  lnterpolant  of  the 
points  (1.1)  of  degree  2n  -  1 . 

II.  If  f(x)  e  Cn_1(I)  is  such  that 

(1.8)  ffx^  -  yL,  (i  -  1,...»k)  , 

(1.9)  f<n-’>(x)  la  absolutely  continuous,  f*n^(x)  8  L2(I)  , 
then 

i  <  .  n  \  (  /  a  (  **)  f  .,1  1 *  f  I 


(1.10) 


/  (f(n,(x))2dx  >  /  (S'n,(x))2dx 
I  I 


with  the  equality  sign  only  if  f(x)  «  S(x)  in,  I. 


'  3 

V;<tv  ^  ^  ->  1  ^ ^wwgp Bia^wwwPBWWWM^MWWIWMi 


III.  If  f(x)  e  d1-^!)  aatigfies  (1.8)  and  (1.9),  than 

(1.11)  /  (f(n))2dx  -  /  (S(n))2dx  -  /  (f(n)  -  S{n))2dx  . 
Ill 

(See  for  instance  (1,  110-116]). 

Our  main  result  is 

Theorem^*  1»  There  is  a  unique  polynomial  Pa(x)  satisfying  (1.6)/ 
where  n  is  defined  by  (1.5). 

2.  We  have 

(1.12)  lim  Pm(  x)  =■  S(x)  uniformly  in  x  0  I  . 

rax*> 

In  view  of  the  extremum  property  II  of  S(x),  the  limit  relation  (1.12) 
may  not  seem  surprising.  Even  so  it  is  no  immediate  consequence  and  our  proof 
of  Theorem  1  occupies  the  remaining  three  sections  of  this  note. 

2.  The  existence  and  uniqueness  of  P(x).  1.  Without  loss  of 

generality  we  may  restrict  the  search  for  P„(x)  to  the  subclass  IT*  c  II  of 

m  mm 

m 

polynomials  P(x)  =*  £  arxr/r!  satisfying 

0 

(2.1)  /  (P(n)(x))2dx  <  +  1  , 

where 

m 

(2.2)  P(n)(x)  =>  l  arxr_n/(r  -  n)  !  . 

r=n 

Let  X^(x),  (i  ■  0,1,...)  be  the  orthonormal  polynomials  for  the  interval  I, 
and  let 


-3- 


(2.3) 


P(n) (x) 


*-n 


I 


CjX^x)  . 


ana 

From  (2.1)  and  FarMTil'i  theorem  we  conclude  that  J  K  » 

0  * 

and  hence  that  |CjJ  <  K,  ( i  -  0,...,m  -  n)  .  Fran  (2.3)  it  followe  that  the 
coefficients 


(2.4)  *n'*a+1' ***'*«  are  bounded  . 

Because  n  <  k  we  have  that  P(x^)  “  y^  for  1  “  1,...,n.  Solving  this 
systesi  for  the  unknowns  a^ , . . .  ,an_ ,  in  terms  of  the  coefficients  (2.4),  we 
conclude  that  for  an  appropriate  constant  I  we  have 

laj  <  H,  (i  -  0,  ...,*)  . 

Now  familiar  compactness  arguments  will  insure  the  existence  of  PM(x) 
satisfying  (1.6)  and  (1.5). 

2.  Let  Pq ( x)  and  p^(x)  be  two  polynomials  in  such  that 

(2.5)  J  (p£n))2dx  -  J  (p\n))2dx  -  . 


Bvidently  also 

(2.6)  pt(x)  -  (1  -  t)p„(x)  +  tp,(x)  •  n^,  (0  <  t  <  1) 

and 

(2.7)  *(t)  -  /  ((1  -  t)pjn)(x)  +  tpjn)(x))2dx  - 


is  a  quadratic  polynomial  in  t  which  is  seen  to  satisfy  the  equations 

(2.8)  *(0)  -  0,  *(1)  -  0  . 

Moreover,  by  (2.7) 

(2.9)  *>(t)  -  t 2  /  (p{n)  -  p^n))2dt  +  At  +  B  . 

I 


Let  us  show  that  the  Inequality 


/  <Pin)  -  Pon))2*t  >  0 


is  impossible.  Indeed,  it  would  imply  by  (2.8)  and  (2.9),  that  l(tQ)  <  0 
for  some  tQ  with  0  <  t,8  <  1.  But  then,  by  (2.7),  we  would  have 
/  ( ^  )2dx  <  n,  contradicting  the  definition  of  m  as  the  minimum. 

j  0  1  9 

We  must  therefore  have 

/  (pjn)  -  pjn))2dx  -  0,  hence  p5n)(x)  «  pjn)(x)  . 


But  then  p^(x)  -  pQ(x)  +  R(x),  where  R(x)  0  ir^^.  Since  R(xi)  •  0  for 
i  “  1,...,k,  and  k  exceeds  the  degree  n  -  1  of  R(x) ,  we  conclude  that 
R(x)  =  0  and  therefore  Pq(x)  *  p^(x). 


3.  Proof  that  pj^(x)  ♦S^n^(x),  as  m  +  »,  in  the  L_-norm.  Let  us  show 


(3.1) 


lim  /  |P<,n)(x)  -  S^(x)|2dx  *  0  . 
mx«  I 


From  the  Property  III,  in  particular  (1.11)  applied  to  f(x)  *  P  (x),  we 


obtain 


(3.2)  /  (P^n)  (x))2dx  -  /  (S'n' (x))2dx  =  /  (P^n)  (x)  -  Stn) (x))2dx 

III 


¥he  definition  (1.5)  of 


(3.3) 


V«"  /  (Pin)(x))2dx 

I 


as  a  minimum,  and  (3.2),  show  that 

(3.4)  /  (pj,n>  (x)  -  S(n>  (x))2dx  -  min  /  (P(n) (x)  -  S(n)  (x))2dx  . 

i  pen  I 

ra 


Clearly,  the  class  11^  expands  on  increasing  m;  this  shows  that  m  is 


-5- 


1 

« 

41 


non-increasing  for  increasing  m,  and  by  (3.2)  also  the  right  side  of  (3.2) 
form  a  non- increasing  sequence.  This  insures  the  existence  of  the  non¬ 


negative  limit 


(3.5) 


lim  /  (P^n)  -  S(n))2dx  -  L  . 
m-H»  I 


A  proof  of  (3.1)  is  now  equivalent  to  showing  that 

(3.6)  L  ■  0  . 

This  requires  two  lemmas  from  Approximation  Theory,  the  first  of  which  is 
well  known  as  an  easy  consequence  of  Weierstrass*  theorem. 

Lemma^l.  Given  e  >  0  we  can  find  a  polynomial  P*(x)  such  that 

(3.7)  |S(x)  -  P„(x)  |  <  e  and  |S(n)(x>  -  pjn)(x)|  <  e  in  I  . 

Indeed,  if  in  the  relation 


n-1  .  .  1  x 

S(x)  *  l  S(r'(a)(x  -  a)  r/r!  +  --  _  }  ^  /  (x  -  t)n"1S(n)  (t)dt 


We  approximate  to  s'  ' ( t)  closely  by  a  polynomial  p(t),  then  the 
polynomial 


n-i  -  x 

p(x)  =  l  S (r)(a)(x  -  a)r/rl  +  --------  /  (x  -  t)n-'p(t)dt 


will  also  approximate  closely  to  S(x).  Since  p'n,(t)  -  p(t),  the  lemma 


follows. 


Lemma  2.  Given  6  >  0,  we  can  find  an  m  >  k,  and  a  P(x)  e  ]T_,  such 

33S3333  ■  1  '  *  III 


(3.8) 


[n)(x)  -  P(n)(x)|  <  5  for  x  e  I  . 


Notice  that  P(x)  e  II  requires  that  P(x,)  -  y* .  This  we  derive  from 

m  1  1 

Lemma  1  by  Lagrange  interpolation  as  follows.  Let  P»(x)  be  the  polynomial 
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- 


of  Lemma  1  satisfying  (3.7),  and  let 


(3.9) 


m  -  maxOt,  degree  of  P„)  , 


hence  P#(x)  e  ir  and  ra  >  k.  Prom  S(x^)  ■  and  the  first  inequality 


(3.7)  we  have 


(3.10)  Ip#^)  -  y^ I  <  e,  (i  -  1, . . . ,k)  , 

Let  Q(x)  e  be  such  that 

(3.11)  Q(xt)  -  P*(x1)  -  yt,  (i  -  1,...,k) 


Finally,  we  define 


(3.12) 


P(x)  =*  P#(x)  -  Q(x)  . 


Notice  that  S(n)(x)  -  P(n)(x)  -  S(n)(x)  -  pin)(x)  +  Qtn)(x)  and  therefore 
(3.13)  |S(n)(x)  -  P(n)(x)|  <  |S(n)(x)  -  Pin)(x)|  +  |Q(n)(x)|  (x  e  I)  . 

At  this  point  observe  that  by  (3.11)  we  have  by  Lagrange's  formula 


Q(x)  *  l  ^(xX P*^)  -  y^ 


and  therefore  also  the  inequality 

(3.14)  maxlP^x^  -  y.^  |  where  *  max  £  |t|n^(x)| 

i  x  i 


By  (3.7),  (3.13),  and  (3.14),  we  conclude  that 


(3.15) 


|S(n)(x)  -  P<n)(x)|  <  (1  +  K  )e  . 

™  n 


Clearly,  P(x)  satisfies  (3.8)  if  we  choose  c  ■  5/(1  +  K  ). 

n 

The  P(x)  defined  by  (3.12)  satisfies  all  conditions  required  by  Lemma 

2s  P(x)  9  n  by  (3.12).  Also  P(x)  &  II  ,  because  by  (3.11) 
m  m 

P(x^)  =  P*(x^)  -  Q(x^)  =  P*(xi)  -  p*(xt)  +  y^  *  y^.  This  completes  a  proof  of 
Lemma  2. 
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A  proof  of  (3*6)  follows  at  once,  because  by  (3.4)  we  have 

0  <  L  <  /  (P^n)  -  S(n))2dx  <  /  (P(n)  -  S(n))2dx  <  62(b  -  a) 

"  I  I 

where  6  is  arbitrarily  small. 

4.  A  proof  of  the  limit  relation  (1.12).  Newton's  formula  with  divided 
differences 

S(x)  =  S(x.,)  +  (x  -  x^JStx^Xj)  +  •••  +  (x  -  x.,)  •••  (x  -  xn_1)S(x1,...,xn) 

+  (x  -  Xl)  •••  (x  -  xn)S(x1,...,xn,x) 

shows  the  following:  If  Qq(x)  e  if  denotes  the  Lagrange  interpolation 
of  S(x)  at  the  points  x^,...,xn,  then 

(4.1)  S(x)  =  Q0(x)  +  (x  -  x.,)  •••  (x  -  xn)S(x1,...,xn,x)  . 

This  is  possible  because  of  (1.2),  hence  n  <  k. 

Now  we  use  the  expression  of  divided  differences  in  terms  of  B-splines: 

If 

(4.2)  M(t)  =  M(t;x1 ,x2, • • ./Xn,x)  (x  e  I) 

is  the  B-spline  of  degree  n  -  1  based  on  the  n  +  1  knots  x1f...,xn,x, 
then 

(4.3)  S(x.,...,x.x)  *  /  M(t)S(n,(t)dt  . 

1  11  n  i  j 

(See  e.g.  [1,  p.  112].  In  that  paper  B-splines  are  still  called  fundamental 
splines.)  Applying  (4.1),  (4.2),  and  (4.3)  to  S(x),  as  well  as  Pm(x),  and 
subtracting  one  equation  from  the  other,  we  obtain  that 

n 

(4.4)  P_(x)  -  S(x)  »  ~  |  |  (x  -  x.)  •  |  M(t)(P^,n)(t)  -  S(n)(t))dt  . 

m  j  j  m 

Applying  the  Schwarz  inequality  we  obtain 


(4.5) 


|P  ( x>  -  S(x)|2  <  (nl)-2  |  |  (x  -  xj2  /  M(t»x1# . .  .,xn,x)2dt  * 

1  J  I 

*  /  (P„n)(t)  -  S(n)(t))2dt  . 

Z 

Since 

— n o  r  9 

(nl)  £  \  I  (x  -  x  )  /  M(t ;x1 , . . . ,xR,x)xdt 

1  3  I 

is  certainly  a  continuous  function  of  the  variable  x  6  I,  it  is  also 

2 

bounded.  Therefore  there  is  a  constant  H  such  that  (4.5)  gives  the 
estimate 

| P  ( x )  -  S(x) | 2  <  H 2  /  (pj,n)(t)  -  S(n)(t))2dt  for  x  e  I  . 

m  x  m 

i 

Now  the  relation  (3.1)  completes  our  proof  of  (1.12). 

5.  Numer ical_exampleg .  The  explicit  evaluation  of  the  polynomial 
Pm(x)  is  an  elementary  problem  of  linear  algebra  in  m  +  1  unknowns.  This 
is  the  reason  why  Theorem  1  is  so  welcome:  It  replaces  for  large  m,  the 
construction  of  Pm(x)  by  the  much  simpler  construction  of  S(x).  We  may  say 
that  Theorem  1  adds  to  the  interest  that  we  attribute  to  the  natural  spline 
interpolant  S(x)  *  s2n-1^* 

The  unicity  of  Pm(x)  in  Theorem  1  clearly  implies  that  if  the  data 
(1.1)  are  symmetric  about  the  origin,  i.e.  b  *  -a,  x^  =  -xk_i+1,  then 
Pm(x)  must  be  an  even  polynomial,  hence  p2r+1^x3  “  P2r(x)* 

For  our  examples  we  choose  the  simplest  such  symmetric  case,  namely 
k  -  3,  (a,b)  =  (-1,1),  x1  =■  -1,  x2  «  0,  x3  »  1,  y1  =  1,  y2  =  0,  y3  =  1  . 
Selecting  n  =  1,  and  m  =  3,  4,  5,  6,  and  7,  we  find  by  elementary 
calculations  that 


9- 


P3(X)  -  X2  , 

P4(X)  -  P5(x)  -  { 18/11 )x2  -  (7/11)x4  , 

Pg(x)  -  P7(X)  -  (25/1 1)x2  -  (25/11 )x4  +  x6  , 
while  the  natural  spline  interpolant  is  the  linear  spline  S(x)  -  |x|, 

-1  <  x  <  1.  The  sequence  of  values 

P3(1/2)  -  .25,  P4(1/2)  -  P5< 1/2)  -  .37,  Pg(1/2)  -  P?(1/2)  -  .44  , 
which  converge  to  S(1/2)  ■  .5,  illustrates  Theorem  1. 
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Let  the  points 

(1)  '  (i  =  lf.rkf  k  >  2)  ,  a  <  x1  <  x2  <  •••  <  x^  <  b  , 

I  =  [a,b]  ,  (-«  <  a  <  b  <  °>) 

be  prescribed.  Furthermore,  let  m  and  n  be  integers  such  that 
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1  <  n  c  k  <  m  , 

and  define  the  polynomial  class 

II  =  (P(x);  P(x)  e  ir  ,  P(x.)  =  y.  ,  (i  =  l,...,k)}  . 
m  mil 

Within  II  we  determine  P  (x)  as  the  solution  of  the  extremum  problem 
mm 


(2) 


/  (P^  (x))2dx  =  minimum  for  P(x)  e  M 
I  m 


Finally,  let  S(x)  =  S_  , (x)  be  the  natural  spline  interpolant  of  degree 
2n  -  1  of  the  k  points  (1) .  Our  main  result  is 

Theorem  1.  1.  There  is  a  unique  polynomial  (x)  which  is  the 

solution  of  the  minimum  problem  (2) . 

2.  We  have 

lim  P  (x)  =  S (x)  uniformly  in  x  e  I  . 

nr*«> 


